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Abstract. This paper studies Monge parameterization, or differential flatness, of control- 
affine systems with four states and two controls. Some of them are known to be flat, and 
this implies admitting a Monge parameterization. Focusing on systems outside this class, 
we describe the only possible structure of such a parameterization for these systems, and 
give a lower bound on the order of this parameterization, if it exists. This lower-bound is 
good enough to recover the known results about "(x, u)-flatness" of these systems, with 
much more elementary techniques, on the order of this parameterization, if it exists. 

Resume. On s'interesse aux parametrisations de Monge, ou a la platitude, des systemes 
affines a quatre etats et deux entrees. Des travaux anterieurs caracterisent ceux de ces 
systemes qui sont "(x, w)-plats", mais on ne sait pas si certains des sytemes restants sont 
plats, ou non. La conjecture est qu'aucun n'est plat, ni Monge-parametrable. Pour ces 
sytemes, on montre que toute parametrisation est d'un type particulier, et on donne une 
borne inferieure sur l'ordre de cette parametrisation, sumsante pour retrouver, de maniere 
beaucoup plus elementaire, le resultat connu sur la "(x, u)-platitude". 



1. Introduction 

In control theory, after a line of research on exact linearization by dynamic state feedback 
[14, 5, 6], the concept of differential flatness was introduced in 1992 in[7] (see also [8, 9]). 
Flatness is equivalent to exact linearization by dynamic state feedback of a special type, 
called "endogenous" [7], but, as pointed out in that reference, it has its own interest, maybe 
more important than linearity. An interpretation and framework for that notion is also 
proposed in [1, 18, 24]; see [16] for a recent review. 

The Monge problem (see the the survey article [25], published in 1932, that mentions 
the prominent contributions [12] and [4], and others) is the one of finding explicit formulas 
giving the "general solution" of an under-determined system of ODEs as functions of some 
arbitrary functions of time and a certain number of their time-derivatives (in fact [25] allows 
to change the independent variable, but we keep it to be time). Let us call such formulas a 
Monge parameterization, its order being the number of time-derivatives. 

The authors of [7] already made the link with the above mentioned work on under- 
determined systems of ODEs dating back from the beginning of 20 th century; for instance, 
they used [12, 4] to obtain, in [22, 17] some results on flatness or linearizability of control 
systems. 
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Let us precise the relation between flatness and Monge parameterizability : flatness is 
existence of some functions — we call this collection of functions a flat output — of the 
state, the controls and a certain number j of time-derivatives of the control, that "invert" 
the formulas of a Monge parameterization, i.e. a solution t t— > (x(t), u(t)) of the control 
system corresponds to only one choice of the arbitrary functions of time appearing in the 
parameterization, given by these functions. Let us call j the order of the flat output. 

Characterizing differential flatness, or dynamic state feedback linearizability is still an 
open problem [10], apart from the case of single-input systems [5, 4]. The main difficulty 
is that the order of a parameterization or a flat output, if there exists any, is not known 
beforehand: for a given system, if one can construct a parameterization, or a flat output, it 
has a definite order, but if, for some integer j, one prove that there is no parameterization 
of order j, then it might admit a parameterization of higher order, and we do not know 
any a priori bound on the possible j's. In the present paper, we consider systems of the 
smallest dimensions for which the answer is not known; we do not really overcome the above 
mentioned "main difficulty", in the sense that we only say that our class of systems does 
not admit a parameterization of order less than some numbers, but the description of the 
parameterization that we give, and the resulting system of PDEs is valid at any order. 

Consider a general control-affine system in R 4 with two controls, where £ G R 4 is the 
state, W\ and w 2 are the two scalar controls and X , X\ and X 2 are three smooth vector 
fields : 

£ = x (z) + w 1 x 1 (z) + w2X 2 (t) . 

In [19], one can find a necessary and sufficient condition on Xq, X\, X 2 for this system 
to admit a flat output depending on the state and control only (j = according to the 
above notations). Systems who do not satisfy this conditions may or may not admit flat 
outputs depending also on some time-derivatives of the control (j > 0) . This is recalled and 
commented in section 2.4 and 5. 

Instead of the above control system, we study a reduced equation (3); let us briefly explain 
why it represents, modulo a possibly dynamic feedback transformation, all the relevant cases. 
Systems for which the iterated Lie brackets of X\ and X 2 do not have maximum rank can be 
treated in a rather simple manner [19, first cases of Theorem 3.1]; if on the contrary iterated 
Lie brackets do have maximum rank, it is well known (Engel normal form for distributions of 
rank 2 in R 4 , see [3]) that, after a nonsingular feedback {w t = /3 l ' (0+ P' A (O w i+ P^ 2 (0 W 2, 
i = 1, 2, with /3 1 ' 1 /? 2,2 — f] 1 ' 2 ^ 2 ' 1 ^ 0), there are coordinates such that the system reads 

£l=Wl, 6 = 7(£l.£2,£3,&) +&Wi , t3=S{€l,&,b,U)+Uwi, ii = W2 (l) 

with some smooth functions 7 and 5. One can eliminate w\ and w 2 and, renaming £1, £2, £3, £4 
as x, y, z, w, obtain the two following relations between these four functions of time : 

y = j(x,y, z, w) + zx , i = 5(x, y, z, w) + wx (2) 

(this can also be seen as a control system with state (x,y, z) and controls w and x). If 
7 does not depend on w, this system is always parameterizable, and even flat (see [19] or 
Example 2.5 below). If, on the contrary, 7 does depend on its last argument, one can, around 
a point where the partial derivative is nonzero, invert 7 with respect to w, i.e. transform 
the first equation into w = g(x, y, z, y— zx) for some function g, and obtain, substituting into 
the last equation, a single differential relation between x,y, z written as (3) in next section. 

Note that (3) also represents the general (non-affme) systems in R 3 with two controls that 
satisfy the necessary condition given in [22, 23], i.e. they are "ruled"; we do not develop this 
here, see [2] or a future publication. 

The paper technically focuses on Monge parameterizations of (3) . The problem is unsolved 
if g and h are such that system (1) does not satisfy the above mentioned necessary and 
sufficient condition. We do not give a complete solution, but our results are more general 
than — and imply — these of [19]. The techniques used in the present paper, derived from 
the original proof of non-parameterizability of some special systems in [12] (see also [21]), 
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are much simpler and elementary that these of [19]: recovering the results from that paper 
in this way has some interest in itself. 

2. Problem statement 

2.1. The systems under consideration. This paper studies the solutions 1 1— ► (x(t),y(t), 
z(t)) of the scalar differential equation 

i = h(x,y,z,X) + g(x,y,z,X)x with A = y — zx (3) 

where g and h are two real analytic functions Ct — > R, Ct being an open connected subset of 
M 4 . We assume that g does depend on A; more precisely, associating to g a map G : £1 — > M 4 
defined by G(x, y, z, A) = (x, y, z, g(x, y, z, A)), and denoting by g± the partial derivative of g 
with respect to its fourth argument, 

<?4 does not vanish on ft and G defines a diffeomorphism fl — ► G(Ct) . (4) 
We denote by O the open connected subset of M 5 defined from SI by : 

(x, y, z, x, y) e n (x,y, z,y — zx) E il . (5) 

From g and h one may define 7 and 5, two real analytic functions G(Cl) — > R, such that 

G _1 (a;, y, z, w) = (x, y, z, j(x, y, z, w)) and S = ho G _1 , i.e. 

w = g(x,y,z,X)^X = ^(x,y,z,w) , (6) 
Ma, 2/, A) = 5(a;, y, z, g(x, y, z, A)) , 5(x, y, z, w) = h(x, y, z, j(x, y, z, w)) . (7) 

Then, one may associate to (3) the control-affine system (1) in R 4 with two controls, that 
can also be written as (2); our interest however focuses on system (3) defined by g and h as 
above. Let us set some conventions : 

The functions 7 and <5: when using the notations 7 and S, it is not assumed that 
they are related to g and h by (6) and (7), unless this is explicitly stated. 

Notations for the derivatives: We denote partial derivatives by subscript indexes. 
For functions of many variables, like <p(u, . . . , M('),t),...,))W) in (10), we use the 
name of the variable as a subscript : p xu (k-i) means d 2 p/ dxdu^ 1 ^ , ip v (t) means 
dip/dv'- 1 ' in (16-b). Since the arguments of g, h, 7, S and a few other functions will 
sometimes be intricate functions of other variables, we use numeric subscripts for 
their partial derivatives : /12 stands for dh/dy, or 34,4,4 for d 3 g/dX 3 . To avoid confu- 
sions, we will not use numeric subscripts for other purposes than partial derivatives, 
except the subscript 0, as in (xo,yo, Zo,Xo,yo) for a reference point. 

The dot denotes, as usual, derivative with respect to time, and ^ the j th time- 
derivatives. 

The following elementary lemma — we do write it for the argument is used repeatedly 
throughout the paper — states that no differential equation independent from (3) can be 
satisfied identically by all solutions of (3) : 

Lemma 2.1. For M e N, let W be an open subset of R 3+2M and R : W — » R a smooth 
function. If any solution (x(.),y(.), z(.)) of system (3), defined on some time-interval 
I and such that (z(t), x(t), . . . , x^ M ^ (i), y(t), . . . y^ M ^(t)) is in W for all t in I, satisfies 
R(z(t),y(t), . . . ,y( M \t),x(t), . . . ,x^ M \t)) = identically on I, then R is identically zero 
on W. 

Proof. For any X e W there is a germ of solution of (3) such that (z(0),x(0), . . . ,x^ M ^(0), 
i/(0), . . . , y( M \0)) = X. Indeed, take e.g. for x(.) and y(.) the polynomials in t of degree M 
that have these derivatives at time zero; Cauchy-Lipschitz theorem then yields a (unique) 
z(.) solution of (3) with the prescribed z(0). □ 
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2.2. The notion of parameterization. In order to give rigorous definitions without taking 
care of time-intervals of definition of the solutions, we consider germs of solutions at time 0, 
instead of solutions themselves. For O an open subset of R™, the notation Cg°(R, O) stands 
for the set of germs at t = of smooth functions of one variable with values in O, see e.g. 
[11]. 

Let k, £, L be some non negative integers, U an open subset of R fc+ ^ +2 and V an open 
subset of R 2i + 3 . We denote byWc Cg°(R,R 2 ) (resp. V C Cg°(R,R 3 ) ) the set of germs of 
smooth functions t ^ (u(t),v(t)) (resp. t ^ (x(t),y(t), z(t)) ) such that their jets at t = 
to the order precised below are in U (resp. in V) : 

U = {(u, v) £ C °°(R, R 2 )|( U (0)X0), . . . , U W (0), u(0), . . . , « (£) (0)) e U}, (8) 

V = {{x,y, z) £ Cq°(R, R 3 ) | (x(0), y(0),z(0), x(0), y(0), . . . ,x^ (0), y^ (0)) € V}. (9) 

These are open sets for the Whitney C°° topology [11, p. 42]. 

Definition 2.2 (Monge parameterization). Let k,£,L be non negative integers, L > 0, 

k < £, and A" = (x , yo ; ^o, £o, Vo, ■ ■ ■ , J/c^) a P om t in f2 x R 2L ~ 2 (fi is defined in (5)). 
A parameterization of order (k, £) at X for system (3) is defined by 



• a neighborhood V of X in O x R 2i 2 , 



• an open subset U C R k+e+2 and 

• three real analytic functions U — ► R, denoted V, Xj 

such that, with U and V defined from U and according to (8)-(9), and 

r : U — > C^°(R, R 3 ) the map that assigns to (u, ti)6W the germ T(m, u) at t = of 

/ x(t) \ ( p(u(t),u(t),---M k) (t),v(t),v(t),---M i) (t)) \ 
t~ y(t) = ^(u(t),«(t),..., «(*)(*),«(*),*(*),..., , (10) 

V *(*) / V x(u(t),u(t),...,uW(t),v(t),Ht),---M e) (t)) J 
the following three properties hold : 

(1) for all (u, v) belonging to U, T(u, v) is a solution of system (3), 

(2) the map T is open and T(U) D V, 

(3) the two maps U — > R 3 defined by the triples (<p u (k) , , Xtt( fc ) ) an d 
(y„c«) , V'wW j X«M ) are identically zero on no open subset of U . 

Remark 2.3 (On ordering the pairs (k,£)). Since u and v play a symmetric role, they can 
always be exchanged, and there is no lack of generality in assuming k < £. This convention 
is useful only when giving bounds on (k,£). For instance, k > 2 means that both integers 
are no smaller than 2. 

Example 2.4. Consider the equation z = y + (y — zx)x , i.e. (3) with g = A, h = y (and 
fl = R 3 ). At any (x ,y , z ,x ,yo,x ,yo) such that Xo + x$ ^ 1, a parameterization of order 
(1,2) is given by : 

v u + ii (1 — v)u + vu 
x = v, y=—— -, z =——— _ (11) 

V + v A — 1 v + v A — 1 

It is easy to check that (x, y, z) given by these formulas does satisfy the equation, point 2 
is true because the above formulas can be "inverted" byu = — z + yx, v = x (this gives the 
"flat output" see section 7), point 3 is true because ipu, ipi), Xu and \v are nonzero rational 
functions. Here, L = 2 and V can be taken the whole set of (x, y, z, x, y, x, y) £ R 7 such that 
x + x 3 7^ 1 and U the whole set of (u, ii, v, v, v) £ R 5 such that v + v 3 ^ 1. 

Example 2.5. Suppose that the function 7 in (2) depends on x, y, z only (this is treated in 
[19, case 6 in Theorem 3.1]). For such systems, eliminating w does not lead to (3), but to 
the simpler relation y — zx = -f(x, y, z). One can easily adapt the above definition replacing 
(3) by this relation. This system y — zx = 7(2;, y, z) admits a parameterization of order (1,1) 
at any (x , y , z , x , yo) such that x + 73(2:0, Vo, zq) 7^0. 

Proof. In a neighborhood of such a point,, the map (x,x,y, z) 1— » (x, x, y, 7(2;, y, z) + zx) is 



MONGE PARAMETERIZATIONS WITH 3 STATE AND 2 CONTROLS 



5 



a local diffeomorphism, whose inverse can be written as (x,x,y,y) i— ► (x,x,y, x(x,x,y,y)), 
thus denning a map \- Then x = u,y = v,z = x( u iU,v,v) defines a parameterization of 
order (1,1) in a neighborhood of these points. 

Remark 2.6. The integer L characterizes the number of derivatives needed to describe the 
open set where the parameterization is valid. For instance, in Examples 2.4 and 2.5, L must 
be taken no smaller than 2 and 1 respectively. Obviously, a parameterization of order (k, £) at 
(x ,yo,Zo,x ,yo, ■ ■ ■ ,x[ L) , y[ L) ) is also, for V > L and any (x L+1 \y ( Q L+1 \ ...,x ( Q L \y ( Q L ] ), 
a parameterization of the same order at (xq, yo, zq, Xo, yo, ■ ■ ■ , x L ^ , y^ L '). 

The above definition is local around some jet of solutions of (3). In general, the idea 
of a global parameterization, meaning that T would be defined globally, is not realistic; it 
is not realistic either to require that there exists a parameterization around all jets (this 
would be "everywhere local" rather than "global") : the systems in example 2.5 admit a 
local parameterization around "almost every" jets, meaning jets outside the zeroes of a real 
analytic function (namely jets such that x + j 3 (x,y,z) ^ 0). We shall not define more 
precisely the notion of "almost everywhere local" parameterizability, but rather the following 
(sloppier) one. 

Definition 2.7. We say that system (3) admits a parameterization of order (k,l) somewhere 
in if there exist an integer L and at least one jet (xq, yo, zq, xq, yo, . . ., :Eq L \ j/q^) G fl x 
M 2L_2 with a parameterization of order (k,£) at this jet in the sense of Definition 2.2 

In a colloquial way this is a "somewhere local" property. Using real analyticity, it should 
imply "almost everywhere local", but we do not investigate this. 

2.3. The functions S, T and J. Given g, h, let us define three functions S, T and J, to 
be used to discriminate different cases. They were already more or less present in [19]. The 
most compact way is as follows : let u>, u 1 and r\ be the following differential forms in the 
variables x,y,z,X : 

0J 1 =dy-zdx, uj = -2g 4 2 dx + (g 4A h 4 - g 4 h 4 . 4 ) uj 1 - g 4A (dz - gdx) , 

T] = dz — gdx — h i uj 1 . ' 

From (4), lo huj 1 t\r\ — 2g 4 2 dx Ady Adz ^ 0. Decompose do>Ao> on the basis uj, lo 1 , ij, dX, 
thus defining the functions S, T and J (we say more on their expression and meaning in 
section 4) : 

{ST \ 
dwAw = - — dAA?7+ — dA A w 1 + J w 1 A r? A w . (13) 
V 2 S4 2 / 

Example 2.8. Le us illustrate the computation of S, T and J on the following three partic- 
ular cases of (3). For each of them, the table below gives the differential forms to and r\, the 
decomposition of dui A u> on u 1 , u, r\, dX and the resulting S, T, J according to (13). System 
(a) was already studied in Example 2.4. 

(a): z = y+(y-zx)x , (b): z = y+(y- zx)(y- (z-l)x) , (c): i = y+ (y-zxfx . (14) 



system 
(14) 


g(x,y,z,X) 


h(x, y, z, A) 


-w/2 
V 


dw Aw 


S,T,J 


(a) 


X 


y 


dx 

dz — Xdx 





0, 0, 


(b) 


X 


y + X 2 


dy — (z — l)dx 
dz — Xdx — 2 Aw 1 


U) 1 A r] A uj 


0, 0, -1 


(c) 


X 2 


y 


dz + 3X 2 dx 
dz — A 2 dec 


\dXNriNuj 


-12, 0, 
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2.4. Contributions and organization of the paper. If S = T = J = 0, i.e. dwAw = 0, 
system (3) admits a parameterization of order (1,2), at all points except some singularities. 
This is stated further as Theorem 4.3, but was already contained in [19]. We conjecture that 
these systems are the only parameterizable ones of these dimensions, i.e. system (3) admits 
no parameterization of any order if (S, T, J) ^ (0, 0, 0), i.e. if du> A lu ^ 0. 

This is unfortunately still a conjecture, but we give the following results, valid if (S, T, J) ^ 
(0,0,0) (recall that k < £, see Remark 2.3) : 

• system (3) admits no parameterization of order (k,£) with k<2ork — £ = 3 
(Theorem 5.4), 

• a parameterization of order (k, £) must come from a solution of the system of PDEs 
£lj (Theorem 5.1), 

• since a solution of this system of PDEs is also sufficient to construct a parameteri- 
zation (Theorem 3.7), the conjecture can be entirely re-formulated in terms of this 
system of partial differential relations. 

Note that this allows one to recover the results from [19] on (x, ^-flatness 1 . See Remark 5.6 
for details. 

The paper is organized as follows. Section 3 is about the above mentioned partial dif- 
ferential system S^'e ■ Section 4 is devoted to some special constructions for the case where 
S = T = 0, and geometric interpretations. The main results are stated in Section 5, based 
on sufficient conditions obtained in Sections 3 and 4, and necessary conditions stated and 
proved in Section 6. Sections 7 and 8 comment on flatness vs. Monge parameterization and 
then give a conclusion and perspectives. 

3. A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS 

This section can profitably be skipped or overlooked in a first reading; the reader will come 
back when needed to this material that might appear, at first sight, somehow disconnected 
from the thread of the paper. 

It defines S^'e an d its "regular solutions", proves that a regular solution induces a param- 
eterization of order (k, £), and that no regular solution exists unless k > 3 and £ > 4. 

3.1. The equation S^'e i regular solutions. For k and £ some positive integers, we define 
a partial differential system in k + £ + 1 independent variables and one dependent variable, 
i.e. the unknown is one function of k + £+ 1 variables. The dependent variable is denoted by 
p and the independent variables by u, it, . . . , u' fe_1 \ x, v, v, . . . , . Although the names 

of the variables may suggest "time-derivatives", time is not a variable here. 

In R k+l+1 with the independent variables as coordinates, let F be the differential operator 
of order 1 

— 2 i~\ £. 2 i~\ 

Q u (t) ^ Q v (l) ^ > 

i=0 i=0 

where the first sum is zero if k < 1 and the second one is zero if £ < 1. 

Let f2 be an open connected subset of R 4 and 7, 5 two real analytic functions f2 — > R 
such that 74 (partial derivative of 7 with respect to its 4 th argument, see end of section 2.1) 
does not vanish on fl. Consider the system of two partial differential equations and three 
inequations : 

p ui k-i)(F Px - 5(x,p,p x ,p xx )) -p xu ( k -i)(Fp- i{x,p,p x ,p xx )) = 0, (a) 
p u (k-i) p xv (t-i) - p xu (k-i) p v (e-i) — , (b) 

£k'i{ PuV-»*0, (c) (16) 

p v a-D ^ , (d) 

71 +l2Px +73 Pxx +74 Pxxx - S ^ . (e) 



1 The term "dynamic linearizable" in [19] is synonymous to "flat" here. 
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To any p satisfying S^'t' we ass °ci a te two functions a and r, and a vector field E : 

j = Pv<A-D _ -Fp + -f(x,p,p x ,p xx ) E = ^ d d 
p u (k-i) ' P u (*-i) ' 9u( fe_1 ) cM^ 1 ) 

We also introduce the differential operator D (see Remark 3.2 on the additional variables 

x,...,x ( - k+e - 1 ')) : 

a-'+'S^+E'^alff- (18) 

i=0 

Definition 3.1 (Regular solutions of )■ A regular solution of system S^'e ls a re& l 
analytic function p : O — > M, with O a connected open subset of R + +1 , such that the image 
of O by (x,p,p x ,Pxx) is contained in f2, (16-a,b) are identically satisfied on O, the left-hand 
sides of (16-c,d,e) are not identically zero, and, for at least one integer K £ {1, . . . , k + £ — 2}, 

£L> K p 7^ (19) 

(not identically zero, as a function of u, . . . , u^ 0-1 ' , x, v, . . . , w^" 1 ), x, . . . , on O x R^). 
We call it K -regular if if is the smallest such integer, i.e. if ED l p — for alii < K — 1. 

Remark 3.2 (on the additionnal variables x, . . . , ^C^- 1 ) in D). These variables appear in 
the expression (18). Note that D is only applied (recursively) to functions of u, . . . , u^ x \ x, 
v, . . . , v( e ~^ only; hence we view it as a vector field in R k+e+1 with these variables as param- 
eters. In fact, D is only used in ED l p, l<i<k + £— 1. This is a polynomial with respect 
to the variables x, x, . . . , x^ with coefficients depending on it, . . . , u^ 1 ^ ,x,v, . . . , ti^ -1 ) via 
the functions p, 7, S and their partial derivatives. Hence ED' l p — means that all these 
coefficients are zero, i.e. it encodes a collection of differential relations on p, where the 
spurious variables x,x, . . . ,x^ no longer appear. Likewise, ED l p ^ means that one of 
these relations is not satisfied. 

Definition 3.3. We say that system E^f admits a regular (resp. K -regular) solution some- 
where in fl if there exist at least an open connected O C R k+e+1 and a regular (resp. 
if -regular) solution p : O — > R. 

Remark 3.4. It is easily seen that p is solution of S^'e ^ an d om y if there exist u and r 
such that (p, a, r) is a solution of 

Fp + rp u (k-i) = l(x,p,p x ,p xx ) Ep = 0,a x = 0, , 2Q . 

FPx + Tp XtU(k -i) = S(x,p,p x ,p xx ) P u (k-D ^ , t x ^ , a ^ 

Indeed, (16) does imply the above relations with a and r given by (17); in particular, t x ^ 
is equivalent to (e) and a ^ to (d); conversely, eliminating a and r in (20), one recovers 
Note also that, with g and /i related to 7 and 5 by (6) and (7), any solution of the 
above equations and inequations satisfies 

Dp x = h(x,p,p x ,Dp- p x x) + g(x,p,p x ,Dp- p x x)x . (21) 

The following will be used repeatedly in the paper : 

Lemma 3.5. If p is a solution of system S^'e an d 

(1) either it satisfies a relation of the typep x — a(x,p) with a a function of two variables, 

(2) or it satisfies a relation of the type p xx = a(x,p,p x ) with a a function of three 
variables, 

(3) or it satisfies two relations of the type p xxx — a(x,p,p x ,p xx ) and 

Fpxx + T Pxxu( k - 1 '> — i ) {x,P,p x ,p xx ), with ip and a two functions of four variables, 
then it satisfies ED l p = for all i > and hence is not a regular solution of S^'e ■ 
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Proof. Point 1 implies point 2 because differentiating the relation p x — a(x,p) with respect 
to x yields p xx — a x (x,p) + p x a p (x,p). Likewise, point 2 implies point 3 : differentiating 
the relation p XjX = a{x,p,p x ) with respect to x yields p xxx = a x (x,p,p x ) + p x a p (x,p,p x ) + 
PxxCt Px {x,p,p x ) while differentiating it along the vector field F + t d/ du^ k ' x "> and using (20) 
yields 

Fp xx +rp xxu (k-i) = j{x,p,p x ,p xx ) a p (x,p,p x ) + 5(x,p,p x ,p xx ) a Px (x,p,p x ). 

Let us prove that point 3 implies ED l p — ED l p x — ED l p xx = for all i > 0, hence the 
lemma. It is indeed true for i = and the following three relations: Dp = ^f(x,p,p x ,p xx ) + 
xp x , Dp x = 5(x,p,p x ,p xx ) + xp xx , Dp xx = ip(x,p,p x ,p xx ) + xa(x,p,p x ,p xx ), that are 
implied by (18), (20) and the two relations in point 3 allow one to go from i to i + 1 
(ED l x = Ex® = and ED l x = Ex^+V = from the very definition of D and E). □ 

3.2. The relation with Monge parameterizations. Let us now explain how a Monge 
parameterization for system (3) can be deduced from a regular solution p : O — > M of S^'e • 
This may seem anecdotic but it is not, for we shall prove (cf. sections 5 and 6) that all 
Monge parameterizations are of this type, except when g and h are such that dui A uj = 
(see (12)-(13)). 

We saw in Remark 3.4 that (16-e) is equivalent to t x ^ 0; let {uq, . . . ,u Q k 1 \xo 1 vq 1 . . ., 
«q £ ^) 6 O be such that t x (uo, . . . , u k *\ xq, t>o> ■ ■ ■ , vjf ^) ^ 0. Choose any (u Q k \ v^p) G 
K 2 (for instance with = 0) such that 

(fe) i (fe-l) W I (fc-l) tr>r>\ 

u ' - cr(u , . . . ,u ',v ,...,v ) w o =r(u ,...,Uo ,zo,t>o, • • • ,t>o )• (22) 
Then, the implicit function theorem provides a neighborhood V of (uo, • ■ • , ,vo,---, v^) 
in M fe+£+2 and a real analytic map ip : V — > K such that <^(mo, ■ • ■ , > v o, ■ ■ ■ , v q ^) = x a 
and 

t(u, u^Mu ■ ■ ■ v&),v, . . . , v«-») = - a{u, u^^v, v^) „«> (23) 
identically on V. Two other maps V — > R may be defined by 

«<*>,«,... - pK---^ (fe " 1) ^(---)^,---,^" 1) ) ; (24) 

X («,...,« (fc) , « (/) ) = Px (u,...,u^M---),v,...M t - 1) )- (25) 
From these <p, ip and x, one can define a map T as in (10) that is a candidate for a 
parameterization. We prove below that, if p is a regular solution of E^'f, then this T is 
indeed a parameterization, at least away from some singularities. The following lemma 
describes these singularities; it is proved in Appendix A. 

Lemma 3.6. Let O be an open connected subset of R k+e+1 and p : O — > M. be a K-regular 
solution of system E^'f , see (16). Define the map n : O x I** — > M. K+2 by 

I p x {u---vS k ~ 1 \x,v---v( e ~ 1 }) \ 
p(u ■ ■ ■ vS k ~ 1 \x, v ■ ■ ■ u^ -1 )) 
tt(« ■ ■ ■ u^- 1 ) ,x,v- v«-V , x ■ ■ ■ x^ ) = Dp{u ■ ■ ■ U ( fe -D ,x,v- , x) 

\D K p{u ■ ■ ■ ,x,v- u^- 1 ) , x ■ ■ ■ x {K) )) 

(26) 

There exist two non-negative integers io < k and jo < £ such that iq + jo = K + 2 and 
( dn dn dn dn \ 

{ 0u(*-*>) ' " ' ' duik- 1 ) ' dvV-*>) ' " " " ' dvV- 1 ) ) ( ' 



is a nonzero real analytic function on O x 



We can now state precisely the announced sufficient condition. Its interest is discussed in 
Remark 5.5. 
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Theorem 3.7. Let p : O — > M, wii/i O C R fc+£+1 open, be a K-regular solution of system 
S^'e , and in, jo be given by Lemma 3.6. Then, the maps <p,ip,x constructed above define a 
parameterization T of system (3) of order (k,£) (see Definition 2.2) at any jet of solutions 
{xo,yo,Zo,x , . . .,x^ K \y , ■ ■ -,vi K) ) such that, for some u , . . - ,u Q k ~ 1: 

(uo, 



,v , 



(fc-i) 

,u , Xo , Vo , ■ ■ ■ , u 

zo =p x (u ,...,u \v ,...,u 



Vn , X ) , 



Vo' = D ' l P( u (h ■ ■ ■ ,u K 1 -',v ,...,v^ c L ',x Q , . . . 0<i<K, t 

the left-hand sides of (16-c,d,e) are all nonzero at {uq, . . . ,u k x \x{),vo, .. .,v^ and the 
function ED K p and the determinant (27) are nonzero at point (ito, • • ■ , u Q k X \xo, ■ ■ ■ , x K \ 
v ,...,v^~ 1] ) e O x R K . 

Proof. Let us prove that T given by (10), with the maps tp, ip, \ constructed above, satisfies 
the three points of Definition 2.2. Differentiating (23) with respect to and yields 
7"x = I, f v (i)T x — —a, hence the point 3 (cr 7^ from (20)). To prove point 1, let 
u{.),v{.) be arbitrary and x(.),y(.), z(.) be defined by (10). Differentiating (10) with respect 
to time, using relations (24) and (25), taking u^- k \t) from (23), one has 

y(t) = Fp + rp^-v +vW{t)Ep + x(t) z{t) , z(t) = Fp x +Tp XtU v.-i)+vW(t)Ep x +x{t)p xx , 

where F is given by (15) and the argument (u(t) . . . u^ 1 ^ (t), x(t),v(t) . . . t^ -1 ) (i)) for Fp, 
Fp x Ep, Ep x , t, p x u (k~i), p u {k-i) and p xx is omitted. Then, (20) implies, again omit- 
ting the arguments of p xx , one has y(t) = j(x(t),y(t), z(t),p xx ) + z(t)x(t), and z(t) = 
6(x(t),y(t), z(t),p xx ) + p xx x(t). The first equation yields p xx = g(x(t),y(t), z(t),y(t) - 
z(t)x(t)) with g related to 7 by (6), and then the second one yields (3), with h related to 6 
by (7). This proves point 1. The rest of the proof is devoted to point 2. 

Let t 1 — ► (x(t),y(t),z(t)) be a solution of (3). We may consider V(u, v) — (x,y,z) (see 
(10)) as a system of three ordinary differential equations in two unknown functions u,v : 

/'-aK...,/-^,...,/- 1 ')/'-^,..,/- 1 ',^,..,/- 1 ') = 0, (29) 



(28) 



,u (fe 1) ,x,v, 
p(u,...,u ( - k ~ 1 \x,v, 



,(<-!) 



= V, (30) 

Pxtu,...,^*- 1 ),^,...,^- 1 )) = z. (31) 

Differentiating (30) K+l times, substituting from (29), and using the fact that ED l p = 
for i < K (see Definition 3.1), we get 

D*p (u(t), . . . , u^i), v(t), vV-V{t),x{t), *«(*)) = ^(*) . 1 < i < K, (32) 

«W ( t ) ED K p (u(t), (t), v(t), . . . , v^-V ( t ), x(t), . . . , (t)) 

rjK+l 

+ D K+1 p («(t), . . . ,uV°-V(t),v(t), vV-VftMt), x^ K+l Ht)) = (*) • (33) 

Equations (30)-(31)-(32) can be written ( z \ 

y 



n(u, ...... 



,u^ k 1] lX ,v,. 



(34) 



V y [K) J 

with 7r given by (26). From the implicit function theorem, since the determinant (27) is 
nonzero, (30)-(31)-(32) yields u^ la \ . . . ,u < - k ^ 1 \ v^ l ~ io \ . . . as explicit functions of 

u, . . . , u( k ~ l0 ~ x \ v, . . . , t;^ - -? 0-1 ), x, . . . ,x( K \ y,...,y( K ) and z. Let us single out these 
giving the lowest order derivatives : 



u 



(fe-io) =/l( u ,..., u (fe-l- 



-1) 



V, . . 



v (t-j ) = p( Uj . . . 7 u (fc-l-i -l) ) Vj . . . j v ( 



-^- 1 \x,...,x^ K \z,y, 



-J'o-l) 



.,y w ), 



,a;( A '),«,y,...,y( A ')). 



(35) 
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Let us prove that, provided that (x, y, z) is a solution of (3), system (35) is equivalent to 
(29)-(30)-(31), i.e. to T(u, v) = (x, y, z). It is obvious that any t ^ (u(t),v(t),x(t), y{t), z(t)) 
that satisfies (3), (29), (30) and (31) also satisfies (35), because these equations were obtained 
from consequences of those. Conversely, let t t— > (u(t),v(t),x(t),y(t),z(t)) be such that 
(3) and (35) are satisfied; differentiating (35) and substituting each time z from (3) and 
(tiC*-*"), <;(*-*>)) from (35), one obtains 

u (k- l0 +t) = _ _ _ ; u(fe _i_ io _i)^ ; ^ v{e -j -i) jXj _ _ _ j x (K+i) j Zj ^ ; 3/*"+*)), i e N, 
w (<-jo+j) = /-2-J ( M , . . . , w (fe-i-io-i) , v , . . . , t;(«-Jo-i) ? . . . ; 0, y, - . . , j e N. 

(36) 

Now, substitute the values of u {k - l °\ v,V°\ v^-^, v& from (36) into (29), (30) and 
(31); either the obtained relations are identically satisfied, and hence it is true that any 
solution of (3) and (35) also satisfies (29)-(30)-(31), or one obtains at least one relation of 
the form (recall that k < £): 

R(U, . . . , uC*- 1 -*- 1 ), v , . . . , „('-*>-!) , Xi . . . ) x (K+l) iZt ^ yiK+l^ = Q _ 

This relation has been obtained (indirectly) by differentiating and combining (3)-(29)-(30)- 
(31). This is absurd because (29)-(30)-(31)-(32)-(33) are the only independent relations 
of order k,£ obtained by differentiating and combining 2 (29)-(30)-(31) because, on the one 
hand, since D K p ^ 0, differentiating more (33) and (29) will produce higher order differential 
equations in which higher order derivatives cannot be eliminated, and on the other hand, 
differentiating (31) and substituting z from (3), from (29) and y from (32) for i = 1 
yields the trivial = because p is a solution of £2'e> see ^ ne proof of point 1 above. 

We have now established that, for (x,y,z) a solution of (3), T(u,v) = (x,y,z) is equiv- 
alent to (35). Using Cauchy Lipschitz theorem with continuous dependence on the param- 
eters, one can define a continuous map s : V — > U mapping a germ (x,y, z) to the unique 
germ of solution of (35) with fixed initial condition (u, . . . , u( k ~' l0 ~ 1 \ v, . . . , = 
(u ,...,u fe 10 1 \ v , . . . , Vq £ Jo ^). Then s is a continuous right inverse of T, i.e. Tos = Id. 
This proves point 2. □ 

3.3. On (non-) existence of regular solutions of system S^'e • 

Conjecture 3.8. For any real analytic functions 7 and S (with 74 7^ 0), and any integers 
k,£, the partial differential system £2'e ( see (16)) does not admit any regular solution p. 

An equivalent way of stating this conjecture is: "the equations ED' l p — 0, for 1 < i < 
k + £ — 2, are consequences of (16)". Note that "ED l p = 0" in fact encodes several partial 
differential relations on p; see Remark 3.2. If 7 and d are polynomials, this can be easily 
phrased in terms of the differential ideals in the set of polynomials with respect to the 
variables u^ 1 ^ , x, v, . . . , v^ 1 ^ with k + £ + 1 commuting derivatives (all the partial 

derivatives with respect to these variables). 

This is still a conjecture for general integers k and £, but we prove it for "small enough" 
k, £, namely if one of them is smaller than 3 or if k = £ — 3. The following statements assume 
k < £ (see remark 2.3). 



2 In other words, (29)-(30)-(31)-(32)-(33), as a system of ODEs in u and v, is formally integrable (see e.g. 
[3, Chapter IX]). This means, for a systems of ODEs with independent variable t, that no new independent 
equation of the same orders (fc with respect to u and i with respect to v) can be obtained by differentiating 
and combining these equations. It is known [3, Chapter IX] that a sufficient condition is that this is true when 
differentiating only once and the system allows one to express the highest order derivatives as functions of the 
others. Formal integrability also means that, given any initial condition (u(0), . . . ,u (fc) (u), v(0), . . .,/'(0)) 
that satisfies these relations, there is a solution of the system of ODEs with these initial conditions. 
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Proposition 3.9. // system S^'e > with k < t, admits a regular solution, then k > 3, I > 4 
and the determinant 



P u (k-1) P u (*-2) 
P xu (k-1) P xu (k-2) 

is a nonzero real analytic function. 



P u (k-3) 



(37) 



Proof. Straightforward consequence of Lemma 3.5 and the three following lemmas, proved 
in appendix B. □ 



Lemma 3.10. Ifp is a solution of system S^'e and either k = 1 or 



; 



P u (k-1) P u (k-2) 
P xu (k-1) P xu (k-1) 

then around each point such that p u (k-i) ^ 0, there exists a function a of two variables such 
that a relation p x = a(x,p) holds identically on a neighborhood of that point. 



Lemma 3.11. Suppose that p is a solution of £^ 

i>k>2, p uik - n ± , 



7,5 



with 



P u (k-1) P u (k-2) 
P xu (k-1) P xu (k-2) 



?0. 



(38) 



// either k = 2 or the determinant ( 37) is identically zero, then, around any point where the 
two quantities in (38) are nonzero, there exists a function a of three variables such that a 
relation p x .x = a{x,p,p x ) holds identically on a neighborhood of that point. 

Lemma 3.12. Let k = I = 3. For any solution p of £3 '3 , in a neighborhood of any point 
where the determinant ( 37) is nonzero, there exist two functions a and ip of four variables 
such that p xxx = a(x,p,p x ,p xx ) and Fp xx + Tp xxu ( k -i, = ip(x,p,p x ,p xx ) identically on a 
neighborhood of that point. 

4. Remarks on the case where S = T = 0. 

4.1. Geometric meaning of the differential form uj and the condition S — T — 0. 

For (x,y, z) such that the set A = {A e M, (x,y,z,\) G ft} is nonempty, (3) defines, by 
varying A in A and x in K, a surface S in [the tangent space at (x,y,z) to] M 3 . Fixing A 
in A and varying x in R yields a straight line S\ (direction (1, z, g(x, y, z, A))). Obviously, 
S = Uaga^a; S is a ruled surface. For each A <G A, let Pa be the osculating hyperbolic 
paraboloid to S along S\, i.e. the unique 3 such quadric that contains S\ and has a contact 
of order 2 with E at all points of S\ . Its equation is 



(y- 



zx - 



A)( 



x + 



zx ■ 



gx 



h))- 



z-gx-h , hj I- 
34 94 Vy 



ZX ■ 



A) =0 



where we omitted the argument (x,y, z, A) of h and g. With ui, uj 1 , r\ defined in (12) and £ 
the vector with coordinates x,y,z, the above equation reads 

(ui, £_) + (h 44 g 4 - #44/14) A + .944/1 (77, i)-h 



- (V,0-a) 



= , 



<Ji' 9i 

that can in turn be rewritten (w 1 , £) — (cj 3 , £) — a = 0, with w 3 and a some differential 
form and function; u>, lo 3 and a are uniquely defined up to multiplication by a non- vanishing 
function; they encode how the "osculating hyperbolic paraboloid" depends on x,y, z and A. 

We will have to distinguish the case when S and T, whose explicit expressions derive from 
(12) and (13): 

5" = 2 54 34,4,4 - 3g 4 ,4 2 , T = 234/14,4,4 - 334,4/14,4, (39) 



3 General hyperbolic paraboloid: (a n i + a 12 Y + a 13 Z) (a 21 x + a 22 Y + a 2a Z) + a al x + a 32 Y + a iz Z + 
a = 0, where the matrix [a 1J ] is invertible and Y, Z stand for y — zx — A, i — gx — h. It contains S\ if 
and only if a 11 = a 31 = a = 0. Contact at order 2 means a 13 = 0, a 33 = -a 12 a 21 /g^, a 32 = —/14a 33 , 
a 22 = ia 21 (#4/i44 — g44/M)/g4 2 , a 23 = ia 21 g44/S4 2 . Normalization: a 12 = a 21 = 1. 
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are zero. From (13), it means that the Lie derivative of u> along d/dX is co-linear to u>, and 
this is classically equivalent to a decomposition u = ku 2 where k ^ is a function of the 
four variables x,y,x,X but a) 2 is a differential form in the three variables x,y,z, the first 
integrals of d/d\. Then, one can prove that the form Cj 3 = u 3 /k and the function a° = a°/k 
also involve the variables x, y, z only. From lj's expression, one can take for instance k = #44 
or k = gg^4 — 2g 2 (they do not vanish simultaneously because g^ does not vanish). Hence 
S = T = if and only if, for each fixed (x, y, z), the osculating hyperbolic paraboloid P\ in 
fact does not depend on A i.e. the surface £ itself is a hyperbolic paraboloid, its equation 
being 

(u;\i)(L0 2 7 i) + (u; 3 ,t)+d a = , (40) 

where £ is the vector of coordinates x, y, z. This yields the following proposition 4 , where the 
functions a , a 1 , a 2 , 6°, b 1 , c°, c 1 of x, y, z are defined by 

3 

CP- = frdx + a 2 ^ 1 - cMz, u 3 = = b°dx + aV - c°dz, a = %- . (41) 
k k k 

Proposition 4.1. If S and T, given by (12)-(13), or (39), are identically zero on fl, then, 
for any (xq, yo, zq, Xq) in Ct, there exist an open set W C R 3 , an open interval /Cl, with 
(xo, t/O) Zq, Ao) € W x I C Cl, and seven smooth functions W — > K denoted by d°, a , a , b , 
b , c , c 1 such that c° + c x A does not vanish on W x I , c 1 ^ — b 1 ^ does not vanish on W , 
and, for (x,y, z, A) € W x /, x e M and 2€l, equation (3) is equivalent to 

X (b 1 (x,y,z)x + a 2 (x,y,z)X — c\x,y,z)ij + (h a (x,y,z)x + a 1 (x,y,z)X — c°(x,y,z)ij + a°(x,y,z) = 0. 

4.2. A parameterization of order (1,2) if S = T = J = 0. It is known [19] that system 
(3) is (a;,u)-flat (see section 7) if S = T = J = 0. For the sake of completeness, let re-state 
this result in terms of parameterization. We start with the following particular case of (3): 

i = k(x, y, z) x X + a(x, y, z) X + b(x, y, z) x + c(x, y, z) with A = y — zx (42) 

where k does not vanish on the domain where it is defined. Note that Example 2.4 was of 
this type with k = 1, a = b = 0, c = y. For short, define the following vector fields: 

9 ^1 d d , d ~ d d n d 
X = c— , X = — h z— — h b— , X = — — h a— , X = k— . 

oz ox ay oz ay oz oz 

Note that, for h an arbitrary smooth function of x, y and z, X°h, X x h, X 2 h, X 3 h also 
depend on x, y, z only. 

Lemma 4.2. System (42) admits a parameterization of order (1,2) at any {xq, yo, zo, xq, yo, 
xq, yo) such that 

Kx Q + K 2 x 3 + (X 1 K-X 3 b+2aK)x 2 + (X 1 a+X O K-X 3 c-X 2 b + a 2 )x +X a-X 3 c^0. (43) 

Proof. From (43), the two vector fields Y = X 2 +xX 3 and Z= [X°+xX 1 , X 2 +xX 3 }+xX 3 
are linearly independent at point (xq, yo, zq, xq, xq). Let then ft, be a function of (x,y, z,x) 
such that Yh = and Zh ^ 0; its "time-derivative along system (42)", given by h = X°h + 
(X 1 h)x+(Yh)X+(dh/dx)x, does not depend on A: it is a function of (x, y, z, x, x); also, since 
Yh = 0, one has Yh = Zh; finally, Zh ^ implies that dh A dh A dx A dx A dx ^ 0. In turn, 
this implies that (x, y, z, x, x) 1— > (h(x, y, z, x), h(x, y, z, x, x), x, x, x) defines a local diffeomor- 
phism at (xo, yo, zo,io,xo)- Let tp and \ t> e the two functions of five variables such that the 

4 We introduced the osculating hyperbolic paraboloid because it gives some geometric insight on u>, S and 
T, but it is not formally needed: Proposition 4.1 can be stated without it, and proved as follows, based on 
(39) (see also [2]): the general solution of S = is a linear fractional expression g = (6 + b A) / (c + c 1 A) 
where b° , b 1 , c°, c 1 are functions of x,y,z only — this is known, for S/(g4,) 2 is the Schwartzian derivative 
of g with respect to its 4 th argument, but anyway elementary — and 94 7^ translates into b^c 1 — b <f 7^ 0; 
then T = yields h = (a + a 1 A + a 2 A 2 ) / (c° + c 1 A) with a , a 1 , a 2 functions of x, y, z. With such g and h, 
multiplying both sides of (3) by c° + c 1 A yields the equation in Proposition 4.1. 
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inverse of that local diffeomorphism is (u, u, v, v, v) i— ► (v, VK^i m, v, v, v), x(u, u, v, v, v), v, v). 
The parameterization (10) is given by : x = v, y = i[)(u, it, v, v, v), z — \(u, u, v, v,v). □ 

Theorem 4.3. If S = T = J = 0, then system (3) admits a parameterization of order (1, 2) 
at any (xo, yo, zo, xq, yo,%o,yo) € (f2 x R 2 ) \ F, where F C £1 x R 2 is closed with empty 
interior. 

Proof. From Proposition 4.1, (3) and (40) are identical. Since dto 2 A w 2 = (see (13)- 
(41)), there is a local change of coordinates (x,y,z) = P(x,y,z) such that Cj 1 = k'dx and 
u 1 = k"(dy — zdx) with k! ^ 0, k" ^ 0. Hence P transforms (40) into (42), for some n, a, b, c. 
Lemma 4.2 gives <p,ip,X defining a parameterization of order (1,2) for this system. Then 
P -1 o ip, p- 1 o ip, P _1 o x define one for the original system (3), or (40); (O x R 2 ) \ F is the 
inverse image by P of the set defined by (43) . □ 

4.3. A normal form if S = T = and J ^ 0. 

Proposition 4.4. Assume that the functions g and h defining system (3) are such that S 

and T defined by (13) or (39) are identically zero on O, and let (xq, yo, z o, -V)) & ^ be such 

that J(x , yo, zo, A ) ^ 0. 

There exist an open set W C R 3 and an open interval I C R such that (xq, yo, zo, Ao) G 
x J c fl, a smooth diffeomorphism P from W to P(W) C R 3 and six smooth functions 

P(W) — ► R denoted k, a, f3, a, b, c such that, with the change of coordinates (x,y,z) = 

P(x,y,z), system (3) reads 

z = k(x, y, z) (y - a(x, y, z) x) (y - (3(x, y, z) x) + a(x, y,z)x + b(x, y,z)y + c(x, y, z) (44) 
and none of the functions k, a — /3, ay and (3$ vanish on W . 

Proof. From Lemma 4.1, we consider system (40). Let P X ,P 2 be a pair of independent 
first integrals of the vector field c 1 (j^ + z-§^j + b 1 -^; from (41), oj 1 , Cj 2 span the annihi- 

lator of this vector field, and hence are independent linear combinations of dP 1 and dP 2 : 
possibly interchanging P 1 and P 2 or adding one to the other, there exist smooth functions 
k 1 ,k 2 ,f 1 ,f 2 such that Co 1 = k l (dP 2 - /MP 1 ), f 1 - f ^ 0, k l + 0, i = 1,2. Now, take 
for P 3 any function such that dP 1 A dP 2 A dP 3 ^ 0; decomposing a) 3 , we get three smooth 
functions p ,^ 1 ,^ 2 such that a) 3 = p° (— dP 3 +p 1 dP 1 + p 2 dP 2 ), p° ^ 0. The change of 
coordinates P = (P l ,P 2 ,P 3 ) does transform system (40) into (44) with 

K= ^!oP-\ a = foP-\ (3 = foP~\ a=p'oP-\ b = p 2 oP-\ c=^oP- 1 . 

k and a — f3 are nonzero because f 1 — f 2 , k 1 and k 2 are. ay and fa are nonzero because the 
inverse images of aydi AdyAdz and (3%dx Adt/Adi by P are dP 1 AdP 2 Ad/ 1 for i = 1,2, that 
are equal, by construction, to da; 1 A a; 1 /(fc 1 ) 2 and do) 2 A w 2 /(fc 2 ) 2 , which are both nonzero 
(the second one because J ^ 0). □ 

Note that (44) is not in the form (3) unless a — z or (3 — z. This suggests, since ay ^ 
and /?3 ^ 0, the following local changes of coordinates A and B, that both turn (44) to a 
new system of the form (3): 

(x, y, z) i-^ A(x, y, z) = (x, y, a(x, y, z)) and (x, y, z) ^ B(x, y, z) = (x, y, (3{x, y, z)) . 

(45) 

These two systems of the form (3) correspond to two choices ft 1 , g 1 and ft 2 , g 2 instead of the 
original ft, g, and they yield, according to (6) and (7), two possible sets of functions 7 and 
S. These will be used in Theorem 6.5; let us give their explicit expression : 

1 i {x,y,z, W ) = W ~Tu {X,V \ Z) ■ S^n^+n^W'V , i € {1,2} (46) 
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with (these are obtained from each other by interchanging a and (3) : 

to 1 ' = (ai + a a.2 + (» + ba) a 3 ) o A^ 1 , m ' = (k a 3 (a — (3)) o ^4 _1 , 
n 1 ' = a 3 o A^ 1 , n 1,1 = (a 2 + ba 3 ) o A^ 1 , n 1,2 = (na 3 ) o A^ 1 , 

(47) 

m 2 '° = (/3i+/3/3 2 + (a + 6/3)/3 3 )oB- 1 , m 2 ^ 1 = ( K /3 3 (/3 - a)) o B" 1 , 
« 2,0 = 03 o B~ x , n 2 ' 1 = (f3 2 + b (3 3 ) o B- 1 , n 2 ' 2 = (k (3 3 ) o B^ 1 . 

Example 4.5. System (14-b) in Example 2.8 is already as in (44). The above choices are, 
for this system: 

7 1 (x,y, z,w) = w, S 1 (x,y,z,w) = y + w 2 , 7 2 (x,y, z,w) = -w, S 2 (x,y,z,w) = y + w 2 . (48) 



We gather here our main results in a synthetic manner. They rely on precise local results 
from other sections : sufficient (sections 4 and 3.2) or necessary (section 6) conditions for 
parameterizability, results on solutions of the partial differential system £jjf'* (section 3.3) and 
on the relation between flatness and parameterizability (section 7). We are not able to give 
local precise necessary and sufficient conditions at a given point (jet) because singularities 
are not the same for necessary and for sufficient conditions; instead, we use the "somewhere" 
notion as in Definitions 3.3 and 2.7. 

Theorem 5.1. System (3) admits a parameterization of order (k,£) somewhere in Q if and 
only if 

(1 ) either S = T = J = on fl (in this case, one can take (k, t) = (1, 2) ), 

1 ^1 2 ^2 . 

(2) or S = T = on Q and one of the two systems £2 [ or [ with 7 1 , 8 l given by 
(46)-(47), admits a regular solution somewhere in tt. 

(3) or S and T are not both identically zero, and the system S^'e w ^ 7 an d & defined 
from g and h according to (6) and (7) admits a regular solution somewhere in £1. 

Proof. Sufficiency : the parameterization is provided, away from an explicitly described set 
of singularities, by Theorem 4.3 if point 1 holds, and by Theorem 3.7 if one of the two other 
points holds. For necessity, assume that there is a parameterization of order (k, £) at a point 
(x, y, z,x,y, . . . ,x( L \y^) in (CI x M. 2L ~ 2 )\F. From Theorems 6.2 and 6.5, it implies that 
one of the three points holds. □ 

Example 5.2. Consider again systems (a), (b) and (c) in (14). From point 1 of the theorem, 
system (a) admits a parameterization of order (1,2), see also Example 2.4. System (b) is 
concerned by point 2 of the theorem: it has a parameterization of order k, I if and only one 
of the two systems of PDEs 

p u (k-i) (Fp x -p- p xx 2 )) - p xu (k-i) (Fp ± p xx )) = p u (k-i) p xv (e-i) - p xu (k-i) p v a-i) = , 
p u (*-i)^0, 2^-1)^0, p + p xx 2 ± p xxx ^ 



admits a "regular solution". Point 3 of the theorem is relevant to system (c) because S 7^ 0: 
(c) admits a parameterization of order k, £ if and only there is a "regular solution" p to 



If Conjecture 3.8 is true, neither system (b) nor system (c) admits a parameterization of any 
order. 



5. Main results 



(49) 




(50) 



Theorem 5.1 gives a central role to the system of PDEs El\ . It makes Conjecture 3.8 
equivalent to Conjecture 5.3 below. Theorem 5.4 states that the conjecture is true for k,£ 
"small enough". 



MONGE PARAMETERIZATIONS WITH 3 STATE AND 2 CONTROLS 



15 



Conjecture 5.3. If dto A u> (or (S,T, J)) is not identically zero on fl, then system (3) does 
not admit a parameterization of any order at any point (jet of any order). 

Theorem 5.4. If system (3) admits a parameterization of order (k,£), with k <£, at some 
jet, then either S = T = J = or k > 3 and £ > 4. 

Proof. This is a simple consequence of Theorem 5.1 and Proposition 3.9. □ 

Remark 5.5. If our Conjecture 3.8 is correct, the systems S^'ji never have any regular 
solutions, and the sufficiency part of Theorem 5.1 (apart from case 1) is essentially void, 
and so is Theorem 3.7. However, Conjecture 3.8 is still a conjecture, and the interest of the 
sufficient conditions above is to make this conjecture, that only deals with a set of partial 
differential equalities and inequalities, equivalent to Conjecture 5.3 below. For instance, if 
one comes up with a regular solution of some of these systems £^' e , this will yield a new 
class of systems that admit a parameterization. 

Remark 5.6 (on recovering the results of [19]). The main result in that reference can be 
phrased : 

" (1) is (x, w)-dynamic linearizable (i.e. (x, «)-flat) if and only if S = T = 
J = " . 

Sufficiency is elementary in [19]; Theorem 4.3 implies it. The difficult part is to prove 
that S = T = J = is necessary; that proof is very technical in [19]: it relies on some 
simplifications performed via computer algebra. From our Proposition 7.4, (x, u)-flatness 
implies existence of a parameterization of some order (k,£) with k < 3 and £ < 3. Hence 
Theorem 5.1 does imply the above statement. 

6. Necessary conditions 

6.1. The case where S and T are not both zero. The following lemma is needed to 
state the theorem. 

Lemma 6.1. 7/(5, T, J) ^ (0,0,0) and system (3) admits a parameterization [<p,ip,x) °f 
order (k, £) at point (xo, yo, zo, . . . , x^\ y^) G R 2i+3 , then ip u (k> is a nonzero real analytic 
function. 

Proof. Assume a parameterization where tp does not depend on vP*\ Substituting in (3) 
yields 

x = h(ip, ip,x,ip- x<f) + g(<p, ip,x,i>- x&v ■ 

Since (p does not depend on u^ k+1 \ differentiating twice with respect to u( fc+1 ) yields 

X u w = VVfc) (hi + g^) , = ip u t,k) 2 (h 4 A + g4Aip) . 

If ip u (k) was zero, then, from the first relation, XuW would too, and this would contradict 
point 3 in Definition 2.2; hence the second relation implies that h^A + 54,4^ is identically 
zero. From point 2 in the same definition, it implies that all solutions of (3) satisfy the 
relation : h^i(x, y,z,y — zx) + .94,4(2;, y,z,y— zx)x — 0. From Lemma 2.1, this implies that 
ft-4.4 and g4A are the zero function of four variables, and hence S = T = J = 0. This proves 
the lemma. □ 

Theorem 6.2. Assume that either S or T is not identically zero on fl, and that system 
(3) admits a parameterization of order (fc, £) at X = (xo, yo, Zo, Xo, yo, ■ ■ ■ 1 ^Kuq^) € f2 x 
R 2L ~ 2 , with k,£,L some integers and f,ip,x defined on U C R k+e+2 . 

Then k > 1, £ > 1 and, for any point (uq, . . . , u\^\ vq, . . . , v^) e U (not necessarily sent 
to X by the parameterization) such that 

ip u (.k)(u , . . .,u ( k \v , . . .,v^ e) ) ^ 0, 
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there exist a neighborhood O of (uo, . . . , u Q k 1 \ ip(u a ■ ■ ■ v^) , v , . . . , v^ e 1 ^) in M. k+e+1 and 
a regular solution p : O —* R of S^'ji , related to ip, ip, x by (23), (24) and (25), the functions 
7 and S being related to g and h by (6) and (7). 

Remark 6.3. The regular solution p is if -regular for some positive integer K < k + £ — 2. If 
L > K, Theorem 3.7 implies, possibly away from some singular values of (xq, yo, zq, xq, yo, ■ ■ ■, 
XQ K \yQ K ^), that system (3) also admits a parameterization of order (k, £) at (x , yo, z ,xo,yo, 
. . . , x Q K \yQ K ^). See also Remark 2.6. 

Proof. Assume that system (3) admits a parameterization (ip, ip,x) of order (k,£) at (xo, yo, 
zo,xo, yo, ■■ ■ , x q L \ Vo^)- Since <p u (.k) does not vanish, one can apply the inverse function 
theorem to the map 

(u,ii, . . .,u {k \v,v, . . . ,w (£) ) i-> (u,.. . ,u (fe_1) ,<£(u, . . . ,u {k) ,v, . . .,v (l) ),v, . . .,v (l) ) 

and define locally a function r of k + £ + 2 variables such that 

tp(u,ii, . . . ,m (/c) , v, v, . . . , = x r(u, ii, . . . , u (fe_1) , x, v, v, . . . ,v w ) = u (fe) . (51) 

Defining two functions p, q by substitution of in ip, x, the parameterization can be 
re-written implicitly as 

V = P(u, u, ... , u (fe_1) , x, v,v,..., v^), 

z = q(u, ii, ... , u( fc_1 ) , x, v, v, . . . , v^), (52) 

u (k) _ r ^ u ) . . . j u( fc_1 ) ; x, V, V, . . . , V^). 

We now work with this form of the parameterization and u,u, . . . , u^ k ^ x \ x,x,x, . . . v,v, . . ., 
v^\v^ +1 \ . . . instead of u, ii, . . . , u^ k ^ x \u^ k \ ... v,v, . . . , v^\v^ l+1 \ .... In order to sim- 
plify notations, let us agree that, if k = 0, the list u,u,... , u^" 1 ) is empty and any term 
involving the index fc — 1 is zero (same with £ — 1 if £ = 0) . Let us also define V and Q by 

V = Fp+rp^-i) +«\((-D +v {l+1) p v i t) , Q = Fq+rq u ^-i)+v^q v ^-i)+v^ +1 \ vW , (53) 

with F given by (15). V and Q depend on it, u, ... , u^ 1 ^, x, v, v, . . . , v^ +1 ^ but not on x; Fp 
and Fq depend neither on x nor on v^ e+1 \ When substituting (52) in (3), using y = V + xp x 
and z = Q + xq x , one obtains : 

Q + xq x = h(x,p,q,X) + g(x,p,q,X)x with A = V + x(p x - q). (54) 

Differentiating each side three times with respect to x, one obtains : 

q x = (h 4 (x,p,q,X) + g 4 (x,p,q,X)x) (p x - q) +g(x,p,q,X), (55) 

= (h 4 A(x,p, q, X) + g Aii (x,p, q, X)x) (p x - qf + 2g A (x,p, q, X){p x - q), (56) 

= {h iAA {x,p,q,X) + g iAA (x,p,q,X)x) (p x - q) 3 + 3g 4A (x,p,q, X)(p x - q) 2 . (57) 

Combining (56) and (57) to cancel the first term in each equation, one obtains (see S and T 
in (39)) : 

(T(x, p, q, X) + S(x, p, q, X)x^j (p x -q) 2 =0. (58) 

The second factor must be zero because, if T + Sx was identically zero as a function of 
u, . . . , u( k ~ x \ x, v, . . . w^" 1 ), then, by Definition 2.2 (point 2), all solutions (x(t),y(t), z{t)) 
of (3) would satisfy T(x, y, z, y — zx) + xS(x, y, z, y — zx) — identically, and this would imply 
that S and T are identically zero functions of 4 variables, but we supposed the contrary. The 
relation q = p x implies 

X = V = Fp + rp u( *-u + v w p v(e -i) + v {l+1) p vW (59) 

and (55) then yields p xx = g(x,p,p x , A), or, with 7 defined by (6), 

A = j(x,p,p x ,p xx ) . (60) 
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Since neither p nor Fp nor r depend on v^ i+1 \ (59) and (60) yield p v (t) — 0, i.e. p is a 
function of u, . . . , u^^ 1 ^, x, v, . . . , u( £ ~ 1 ) only. Then (59) and (60) imply (109) with / = 7. 
Furthermore, since <p v (e> 7^ (point 3 of Definition 2.2), (51) implies r v (e) 7^ 0. Also, if p was 
a function of x only, then all solutions of (63) should satisfy a relation y(t) = p(x(t)), which 
is absurd from Lemma 2.1. We may then apply lemma C.l and assert that k > 1, I > 1, 
p u (k-i) ^ 0, p„(<-i) 7^ 0. 

Since p does not depend on , (60) implies that the right-hand side of (59) does not 
depend on «W either; since p u (k-n 7^ 0, r must be affine with respect to i.e. 

r = t + (ju (£) , (61) 

with a and r some functions of u, . . . , u^ 0-1 ), ir, w^" 1 ^. Since p, q = p x , A and q x = p xx 

do not depend on v^'\ (54) implies that Q does not depend on «W either; with p x = and 
r given by (61), the expression of Q v (e) is op xu (k-v> +P xv (t~i) while, from (59), the expression 
of V v (£> . Collecting this, one gets 

crp u ( k -i) + p v (t-i) = ap xu(k -i) +p xv (t-i) = . (62) 

Since p u (k-i) 7^ and p v (i-n 7^ 0, (62) implies Eg = 0, and also <r x = 0, a 7^ 0. 
Then, since r x 7^ (see (51)), (61) implies t x 7^ 0. With the above remarks, (59) yields 
V = A = Fp + rp u (k-i) and hence, from (60), the first relation in (20). In a similar way, 
(53) yields Q = Fp x + rp xu (k-i), and substituting in (54), one obtains (the terms involving 
x disappear according to (60)) Fp x — 6(x,p,p x ,p xx ) + p xu <.k-i)T = with S defined by (7). 
This proves that p satisfies (20), equivalent to (16) according to Remark 3.4, and hence that 
p is a solution of S^.'\. 

To prove by contradiction that it is if-regular for some K < k + £ + 1, assume that 
ED l p = for < i < k + I. Then p x ,p,..., D^+^p, x,..., x ( k+e -^ are 2k + 21 + 1 
functions in the 2k + 2£ variables u, . . . , u^^ 1 ^, v, . . . , v^ _1 \ x, . . . , a;(' c +^ _1 ) . At points where 
the Jacobian matrix has constant rank, there is at least one nontrivial relation between them. 
From point 2 of Definition 2.2, this would imply that all solutions of system (3) satisfy this 
relation, say R(z(t), y(i), . . . , y( fe+£_1 ) (t), x{t), . . . , x ^ k+e ~ 1 ^ (t)) = 0, which is absurd from 
Lemma 2.1. □ 

6.2. The case where S and T are zero. Here, the situation is slightly more complicated: 
we also establish that any parameterization "derives from" a solution of the system of PDEs 
(16), but this is correct only if J is not zero, and there are two distinct (non equivalent) 
choices for 7 and 5. If J 7^ 0, we saw, in section 4, that possibly after a change of coordinates, 
system (3) can be written as (44), which we re- write here without the tildes: 

z = n(x,y,z) (y - a(x,y,z)x) {y - (3(x, y, z) x) + a(x,y,z)x + b(x,y,z)y + c(x,y,z) , (63) 

where k, a, (3, a, b, c are real analytic functions of three variables and k 7^ 0, a — (3 7^ 0, 
da/dx 7^ 0, dfi/dx 7^ 0. We state the theorem for this class of systems, because it is simpler 
to describe the two possible choices for 7 and 5 than with (3), knowing that S = T = 0. 

Lemma 6.4. If system (63) admits a parameterization (<p,ip,x) °f order (k,£) at a point, 
then ip u (k) is a nonzero real analytic function. 

Proof. After a change of coordinates (45), use Lemma 6.1. □ 

Theorem 6.5. Let (xq, yo, zq) be a point where k, a — (3, a>3 and Ps are nonzero, and k, £, L 
three integers. If system (63) has a parameterization of order (k,£) at X = {xq, y$, zq, xq, yo, 
. . . ,x^\ y^) with ip,ip,x defined on U C M. k+e+2 , then k > 1, £ > 1 and, for any point 
(uq, . . . , UQ k \v , . . . , v^) G U (not necessarily sent to X by the parameterization) such that 

if u (k){u , . . .,u ( k \v , . . .,v^ e) ) ^ 0, 
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there exist a neighborhood O of (uq, . . . , 1 \ ip(u a ■ ■ ■ v^) , Vq, . . . , Vq ^) in W k+e+1 and 

1 S 1 2 5 2 

a regular solution p : O — > M of one of the two systems £^ £ ' or £^1 with j % , 8 l given by 
(46)-(47), such that p,tp,ip,x are related by (23), (24) and (25). 

Remark 6.3 applies to this theorem in the same way as theorem 6.2. 

Proof. Like in the beginning of the proof of Theorem 6.2, a parameterization (<p, ip, \) of 
order (k,£) with f u {k) ^ yields an implicit form (52). Substituting in (63), one obtains 
an identity between two polynomials in and x. The coefficient of (v^ e+1 ^) 2 in the 

right-hand side must be zero and this yields that p cannot depend on t)W; the linear term 
in ti^ +1 ) then implies that q does not depend on either. To go further, let us define, as 
in the proof of Theorem 6.2, 

V = Fp + rp uik -i } , Q = Fq + rq u(k -i) +v (e) q v (i-i) , (64) 

with F as in (15). Still substituting in (63), the terms of degree 0, 1 and 2 with respect to 
x then yield 

Q = n{x,p, q)V 2 + b{x,p, q)V + c{x,p, q) , 

q x = K(x,p,q) (2p x - a(x,p,q) - (3(x, p, q)) V + a(x,p, q) + b(x,p,q)p x , (65) 
= (p x - a(x,p,q)) (p x - (3(x,p,q)) . 

The factors in the third equation cannot both be zero because a — (3 7^ 0. Let us assume 

p x - a(x,p,q) = , p x - f3(x,p,q) ^ (66) 

(interchange the roles of a and (3 for the other alternative). Since a>3 ^ 0, the map A defined 
in (45) has locally an inverse A -1 , and the equation in (66) is equivalent to (x,p, q) = 
A^ 1 (x 7 p,p x ); by differentiation an expression of q x as a function of x,p,p x ,p xx is obtained; 
solving the second equation in (65) for V and substituting q and q x , one obtains V = 
^ x (x,p,p x ,p xx ) with 7 1 defined by (46)-(47). If one had chosen the other alternative in (66), 
A and 7 1 would be replaced by B and 7 2 . 

Since V is also given by (64), the relation (109) holds with / = 7 ; also, for the same 
reasons as in the proof of Theorem 6.2 (two lines further than (60)), r„(<o is nonzero and it 
would be absurd that p depends on x only. One may then apply Lemma C.l and deduce 
that k > 1, £ > 1, p u (k-i) 7^ 0, p v (t-i) 7^ 0. 

Since neither p nor V = 7 1 (x,p,p x ,p a;a ;) depend on and p u (k-i) 7^ 0, the first equation 
in (64) implies that r assumes the form (61) with a and r some functions of the k + £ + 1 
variables u,u,..., u^ x \x, v, v, . . . , and that two relations hold: on the one hand 

o~p u (k-i) +p v (e-i) = 0, i.e. one of the relations in (20), and on the other hand the first relation 
in (20) with 7 = 7 1 . Similarly, the second equation in (64) yields aq u (k-i) + q v {i-i) = 
and Fq + rq u {k-i) = Q = kV 2 + bV + c. Applying F + rd/du^-^ and E to the first 
relation in (66) and using the four relations we just established, one obtains on the one hand 
the second relation in (20), with S = S 1 (S 1 defined in (46)- (47)) and on the other hand 
cr Pxu( k - 1 > +£W«-!> — 0- The relations a x = 0, a 7^ and t x 7^ are then obtained exactly 
like at the end of the proof of theorem 6.2; hence p satisfies (20) with 7 = 7 1 and 6 = S 1 ; 
this proves, thanks to Remark 3.4, that p is a solution of £^ [ (it would be £^ £ ' if one 
had chosen the other alternative in (66)). The last paragraph of the proof of Theorem 6.2 
can be used to prove that this solution is if- regular with K < k + £ + 1 . □ 

7. Flat outputs and differential flatness 

Definition 7.1 (flatness, endogenous parameterization [7]). A pair A— (a, b) of real analytic 
functions on a neighborhood of (x , yo, zq, . . . , a;^ , y^) infix M 2 ^ 2 is a flat output of order 
j at X = (xQ,yo, zq, . . . , Xq L \ y^) (with L > j > 0) for system (3) if there exists a Monge 
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parameterization (10) of some order (k, £) at X such that any germ (x(.),y(.), z(.), u(.), v(.)) € 
VxW (with U, V possibly smaller than in (10)) satisfies (67) if and only if it satisfies (68): 

v(u(t),u(t),---M k) (t),v(t)Mt),---M e) (t)) = x(t) ) 

iP(u(t),u(t),---,uW(t),v(t),v(t),---M e) (t)) = V(t) \ , (67) 
X {u(t),u(t),...,u^(t),v(t),v(t),...,vW(t)) - z(t) J 

z(t) = h(x(t),y(t),z(t), y(t)-z(t)x(t)) + g(x(t),y{t), z(t), y(t)- z(t)x(t)) x(t) ] 

u(t) = a(x(t),y(t), z{t), x(t),y(t), x(t),y(t), x& (t), y& (t)) \ . (68) 

v(t) = b((x(t), y(t),z(t),x(t),y(t), x(t),y(t), x& (t), yCi) (t)) J 

System (3) is called flat if and only if it admits a flat output of order j for some j € N. A 
Monge parameterization is endogenous 5 if and only if there exists a flat output associated 
to this parameterization as above. 

In control theory, flatness is a better known notion than Monge parameterization. For 
general control systems, it implies existence of a parameterization (obvious in the above 
definition), and people conjecture [10] that the two notions are in fact equivalent, at least 
away from some singular points. In any case, our results are relevant to both: systems (3) 
that are proved to be parameterizable are also flat and our efforts toward proving that the 
other ones are not parameterizable would also prove that they are not flat. 

Theorem 3.7 gave a procedure to derive a parameterization of (3) from a regular solution 
p of £2'e> an ^ we saw m Section 5 that, unless S = T = J = 0, these are the only possible 
parameterizations. One can tell when such a parameterization is endogenous: 

Proposition 7.2. Let p : O — > IK, with O C R fe+ ^ +1 open, be a regular solution of system 
£^'e ■ The parameterization of order (k,£) of system (3) associated to p according to Theo- 
rem 3. 7 is endogenous if and only if p is exactly (k + £ — 2) -regular; then, the associated flat 
output is of order j < k + £ — 2. 

Proof. In the end of the proof of Theorem 3.7, it was established that (67), written T(u, v) — 
(x,y,z), is equivalent to (35) if (x,y, z) is a solution of (3). If either i < k or j < £ in 
(35), then there are, for fixed x(.),y(.), z(.)), infinitely many solutions u(.), v(.)) of (35) while 
there is a unique one for (68). Hence io = k and jo = £ if (67) is equivalent to (68); then 
K = io + jo — 2 = k + £ — 2 so that p is (k + £ — 2)-regular and (35) (where u and v do not 
appear in the right-hand side) is of the form (68) with j — K = k + £ — 2. □ 

The main result in [19] is a necessary condition for "(x, u)-dynamic linearizability" ((x, u)- 
flatness might be more appropriate) of system (1). For system (1), it means existence of a flat 
output whose components are functions of £ 2 , £ 3 , £ 4 , w 1 , w 2 ; for system (3), it translates 
as follows. The functions 7 and S in (1) are supposed to be related to g and h in (3) according 
to (6) and (7). 

Definition 7.3. System (1) is "(x, u)-dynamic linearizable" is and only if system (3) admits 
a flat output of order 2 of a special kind : A(x, y, z, x, y, x, y) = a(x, y, z, A, x, A) for some 
smooth a. 

The following proposition is usefull to recover the main result from [19], see Remark 5.6. 

Proposition 7.4. If system (1) is "(x,u) -dynamic linearizable" in the sense of [19], then 
(3) admits a parameterization of order (k,£) with k < 3 and £ < 3. 

Proof. I a(x,y,z,X,x,X) \_ 

Consider the map (x, y, z, A, x, A, . . . , x^ , A < 4 > ) h-> .. . a(x ' y > !' } Xl A ' x ' ^ f ,. . 

a(x, y, z, A, x, A, . . . , x w , A w ) 

\a( 3 )(x,y,z,A,x,A,...,x( 4 ),A( 4 ))y 



5 This terminology (endogenous vs. exogenous) is borrowed from the authors of [7, 15]; it usually qualifies 
feedbacks rather than parameterizations, but the notion is exactly the same. 
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Its Jacobian is 8 x 12, and has rank 8, but the 8x8 sub-matrix corresponding to derivatives 
with respect to x, A, ... , , A' 4 ) has rank 4 only. Hence x, y, z, and A can be expressed 
as functions of the components of o, d, d, a^ 3 \ yielding a Monge parameterization of order at 
most (3, 3). □ 

8. Conclusion 

Let us discuss both flatness (see Section 7) and Monge parameterization. For convenience, 
assume k < I and call F-systems the systems (3) such that S = T = J = and C-systems 
all the other ones. 

F-systems are flat; this was proved in [19]. This paper adds that they admit a Monge 
parameterization of order (1,2), but does not prove differential flatness of any system not 
known to be flat up to now: C-systems are not believed to be flat. It does not either prove 
non-flatness of any system: it only conjectures that no C-system admits a parameterization, 
and hence none of them is flat. To the best of our knowledge, no one knows whether simple 
systems like (14-b) or (14-c) are flat of not. 

The first contribution of the paper is to prove that a C-system admits a parameterization 
of order (k,£) if and only if the PDEs S^'f , for suitable 7, 5, admit a "regular solution" p. 
The second contribution is to prove that, for any 7, S, there is no regular solution to S^'e ^ 
either fc<2or/c = £ = 3 (this does not contradict existence of parameterizations of order 
(1,2) for F-systems: these do not "derive from" a solution of these PDEs). We guess, in 
Conjecture 3.8, that even for higher values of the integers k,£, none of these PDEs have any 
regular solution; this would imply that C-systems are not flat. 

Besides recovering the results from [19] with far more natural and elementary arguments, 
we believe that some insight was gained on Monge parameterizations of any order for "C- 
systems", by reducing non-parameterizability to non-existence of solutions to a systems of 
PDEs that can easily be written for any k,£. 

The main perspective raised by this paper is to prove Conjecture 3.8. The only the- 
oretical difficulty is, in fact, that no a priori bound on the integers k,£ is known. Indeed, 
as explained in Section 3.3, for fixed k,£, 7, 5, it amounts to a classical problem. To prove 
Proposition 3.9, we solved, in a synthetic manner, that problem for fc<2orfc = i? = 3 
and arbitrary 7 and (5. We lack a non-finite argument, or a better understanding of the 
structure, to go to arbitrary k, I. Let us comment more on the (non trivial) case where 7 
and S are polynomials, for instance the very simple ones in (49). For fixed k,£, the question 
can be formulated in terms of differential polynomial rings: does the differential ideal gen- 
erated by left-hand sides of the equations (49) contain the polynomials ED l p ? Differential 
elimination (see [20] or the recent survey [13]) is relevant here; finite algorithms have been 
already implemented in computer algebra. Although we have not yet succeeded (because of 
complexity) in carrying out these computations, even on example (49) for (k,£) = (3,4), and 
although it will certainly not provide a bound on k, £, we do believe that computer algebra 
is a considerable potential help. 

Another perspective is to enlarge the present approach to higher dimensional control 
systems. For instance, what would play the role of our system of PDEs S^'e when, instead 
of (3) , one considers a single relation between more than three scalar functions of time (this 
captures, instead of (1), control affine systems with n states and 2 controls, n > 4) ? We 
have very little insight on this question: the present paper strongly takes advantage of the 
special structure inherent to our small dimension; the situation could be far more complex. 
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Appendix A 

Proof of Lemma 3.6. For this proof only, the notation !F% t j (0 < i < k, < j < £) stands 
either for the following family of i+j vectors in ~R K+2 or for the corresponding (K+2) x (i+j) 
matrix : 

/ dir dir dir dir 

iJ = \du( k ->) ' ' ' ' ' 9u( fe -!) ' dv«-i) ' ■ ■ ■ ' dv(t-V , 
with the convention that if i or j is zero the corresponding list is empty; T%j depends on 
u, . . . , u( fc_1 ), v, . . ., v^ 1 ), x, . . . , a;^. Let us first prove that, at least outside a closed subset 
of empty interior, 

Rank.F M = K + 2 . (69) 

Indeed, if it is smaller at all points of O x R K , then, around points (they form an open 
dense set) where it is locally constant, there is at least one function R such that a non-trivial 
identity R(p x ,p, . . . , D K p, x, . . . , x^) = holds and the partial derivative of R with respect 
to at least one of its K + 2 first arguments is nonzero. Since p is isT-regular, applying E to 
this relation, shows that R does not depend on D K p, and hence does not depend on x^ 
either. Then, applying ED, ED 2 and so on, and using the fact that, according to (21), 
Dp x is a function of p x ,p, Dp, x, x, we get finally a relation R(p x ,p,x) = with (R P:c ,Rp) ^ 
(0,0). Differentiating with respect to u^ x \ one obtains R Px P xu (k-i) + R p p u (k-i) = 0; 
hence, from the first relation in (16-c), i? Px ^ 0, and the relation R(p x ,p,x) — implies, 
in a neighborhood of almost any point, p x = f(p, x) for some smooth function /. From 
Lemma 3.5, this would contradict the fact that the solution p is X-regular. This proves (69). 

Let now W s (1 < s < K + 2) be the set of pairs such that i + j = s and the rank of 
Ti.j is s at least at one point in O x R K , i.e. one of the s x s minors of !Fij is a nonzero real 
analytic function on O x R^. The lemma states that Wk+2 is nonempty; in order to prove 
it by contradiction, suppose that Wk+2 = and let s be the smallest s such that W s = 0. 
From (16-c), Wi contains (1,0), hence 2 < s < K + 2 < k + l+1. Take (i',f) in Ws-i; Tv ,y 
has rank i' + j' (i.e. is made of i' + j' linearly independent vectors) on an open dense set 
A C O x R K . Let the i\ < k and ji < I be the largest such that Ti x j* and Ty , 3 \ have rank 
s — 1 on A. On the one hand, since i' + j 1 = s — 1 < k + 1, one has either i' < k or j' < I. 
On the other hand since W-g is empty, it contains neither [i 1 + 1, j') nor + 1); hence the 
rank of -Tv+i .y is less than i' + j' + 1 if i' < k, and so is the rank of Ty ,y+\ if j' < £. 

To sum up, the following implications hold: i' < k => i\ > i' + l and j' < £ +> ji > j'+l ■ 
From (69), one has either ii < k or ji < £. Possibly exchanging u and v, assume i\ < k; all 
the vectors dir / du^ k ~ %1 \ . . . , dn/du^- k ~ l +1 \ dn / 'du^ e ~^\ . . . , dir/du^i +1 ) are then linear 
combinations of the vectors in Fi>j>, while dir j du^ k ~ n ~ 1 " > is not : 

RankJT;, j, = i'+f , Rank^ n dl = i'+f , Rank (^-^— - , jr.,^ = i'+j' + l (70) 

on an open dense subset of O x R K , that we still call A although it could be smaller. In 
a neighborhood of any point in this set, one can, from the third relation, apply the inverse 
function theorem and obtain, for an open f2 C ]H k + i + K + 1 ) a map £1 — > W +J +1 that ex- 
presses u( k ~ l . ,u( k ~ x \ \ . . . ,v( £_1 ) and u^ - * 1-1 ) as functions of u, . . . ,u( k ~ ll ~ 2 \ 
u^ k ~ ll \ . . . ,u( k ~ l v, . . . ,v( l ~i ~ x \ x, . . . ,x( K ~) and i' + j' + 1 functions chosen among 
p x ,p, Dp, . . . , D K p (i'+j'+l columns defining an invertible minor in (97r/9zt^ _n_1 ^ , Fi'j')). 
Focusing on one has 
u (fe-»i-i) = 



B 



(u,..., u^ k - il - 2 \u^- il \ u^- 1 '- 1 ^^, i^-j"- 1 ), x, x^ K \p x ,p, D K p 



(71) 

where B is some smooth function of k + £ + 2K + 2 — i' — j' variables and we have written 
all the functions p x ,p, Dp, . . . , D K ~ 1 p although B really depends only on i' + j' + 1 of them. 
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Differentiating (71) with respect to u( k ~ l \ . . . , u^ -1 ) , \ . . . ,v ( - i ~ 1 \ one has, with 
obvious matrix notation, fj^ ■ ■ ■ jHW^p) ^'j' = ^> where the right-hand side is a 

line- vector of dimension i' +f; from (70), this implies 



dB dB dB 



d Px dp dD K -^ 



Fh,h = 0, (72) 



where the right-hand side is a now a bigger line- vector of dimension i\ + j\. Differentiating 
(71) with respect to u {k ~ ll \ . . . , u^"'' -1 ), v^~ jl \ . . .,v( e ~ j '-V and using (72) yields that B 
does not depend on its arguments u^ k ~ ll \ . . . ,u( k ~ l _1 ) and v^~^ \ . . . ,v^~i ~ x \ B cannot 
depend on D K p either because ED K p ^ and all the other arguments of B are constant 
along E; then it cannot depend on x^ either because appears in no other argument; 
(71) becomes 

u (k-i!-i) = B (u, . . .,u^~ 2 \v, . . a, ■ ■ -^x^^Kp^P, • • ■,D K - 1 p) . 

Applying D, using (21) and substituting u( fc_Jl_1 ) from above, one gets, from some smooth 

C, ' 

u< k -^ =C(u,.. .,u^- 2 \v, . . .,v^\x, . . .,x {K \p Xl p, . . .,D K p) . (73) 

v ' 

empty if ji =£ 

Differentiating with respect to u^ k ~ l \ , u^ fe_1 \ \ , i/^" 1 ) yields 
(jjF" T5p ' ' ' dD^-ip ) ^ v <i' = ^' ^ ne r ight-hand side being a line-vector of dimension i' 
From the first two relations in (70), dir / (9u( fc_ll_1 ) is a linear combination of the columns of 

f dC dC_ dC \ dir 

\dp~x~ ~dp~ " ' dD K ~ 1 p J du^-ii-V 
This implies that the derivative of the right-hand side of (73) with respect to u^ - * 1-1 ) is 
zero. This is absurd. □ 

Appendix B. Proof of Lemmas 3.10, 3.11 and 3.12 

We need some notations and preliminaries. With F, E and r defined in (15) and (17), 
define the vector fields 

X = L Y = F + T J^ < 74 > 
X 1 = [A, Y], X 2 = \X U Y\ , E 2 = [E, Y], E 3 = [E 2 , Y] . (75) 

Then (20) obviously implies 

Yp = l(x,P,Px,Px,x), YPx = S(x,p,p x ,p XtX ) , A<r = 0, (76) 
and a simple computation yields (we recall E from (17)) : 

Ai = t x „ , A 2 = t x „ (f_ 2 ) +('••) fl,,(fc-i) ) (77) 



Ty t j> , hence one also has ( — — • • • ^—^ , ) „ - — — = . 



X Q u (k-1) 1 ^-2 — T X Q u (k-2) "TV'/ a„(fe-l) ' 

d i d d I d , / \ d 

j(t-i) ° 9ttC=-i) j ^2 — g v (e-2) t o g u (k-2) > V ' ' ) g u (k-i) > 

= dv (t-3) + a Q U ii-3) + (• • • ) a„<fc-2) + (' ' ' ) au(fc-i) • 



The vector field Ai and X 2 are linearly independent because t x ^ 0, see (20). Computing 
the following brackets and decomposing on Ai and A2, one gets 

[A, Ax] = AAi , [Ax , A 2 ] = A'Ai + A" A 2 , (78) 
[A, E] = , [A, £ 2 ] = A*Ai , [A, £3] = M' *i + m" A 2 , [£7 2 , A 2 ] - + v" X 2 .(79) 
for some functions A, A', A", /i, jj! , fx", u', v" . 
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Proof of Lemma 3.10. From (16-c), y = p{u, . . . , u^ -1 ), x, v, . . . , i^ -1 )) defines local coor- 
dinates u, . . . tS k ~ 2 } , y,x, v, . . . , . Composing p x by the inverse of this change of coordi- 
nates, there is a function a of k+£+l variables such that p x — a(u, . . . , tS k ~ 2 ^ , p , x, . . . , u^ -1 )) 
identically. Since Ep = Ep x = (see (20)), applying E to both sides of this identity yields 
that a does not depend on its argument v^ x \ Similarly, if k > 2, differentiating both sides 
of the same identity with respect to u^ k ~^ and u^ k ~ 2 \ the fact that the determinant in the 
lemma is zero implies that a does not depend on its argument w/ fe ~ 2 ). To sum up, p and p x 
satisfy an identity 

p x = a(u, u (fc ~ 3) ,p,x,v,..., v {l ' 2) ), 
where the first list is empty if k = 1 or k = 2. Now define two integers m < k — 3 and n < I— 2 
as the smallest such that a depends on u, . . . , u^ m \x, y, v, . . . , v^ n \ with the convention that 
m<0iffc = l,fc = 2, or a depends on none of the variables u, . . . , u( k ~ 3 ^ and n < if a 
depends on none of the variables v, . . . , v^ 2 \ 

Applying Y to both sides of the above identity yields 

m n 

Y Px = a v Yp + J2 uit+1)a u(» +Xy i+1)c V<) , 

i=0 i=0 

where, if m < or n < 0, the corresponding sum is empty. Using (76), since p xx — a x + aa y , 
one can replace Yp with 7(21, y, a, a x + aa y ) and Yp x with S(x, y, a, a x + aa y ) in the above 
equation, where all terms except the last one of each non-empty sum therefore depend on 
u, . . . , u( m \ x, y, v, . . . , only. Differentiating with respect to u( m+1 ) and yields 
a u (») = a v (n) — 0, which is possible only if m < and n < 0, hence the lemma. □ 

Proof of Lemma 3.11. From (38), setting 

y = p{u, u^lx, v,..., vV-V), z = Px (u, u^lx, v^) 7 (80) 

one gets some local coordinates (u, . . . , u^ k ~ 3 \ x, y,z,v, . . . , v^ 1 " 1 ). In these coordinates, the 
vector fields X and Y defined by (74) have the following expressions, where X an d ot are 
some functions, to be studied further : 

Odd , . 

X = ^r + zj- + a— , 81 
ox ay Oz 

Y - &&*J*+Z*<» n m+'t** n m- (82) 

In the expression of Y, the third term is zero if k = 2, the fourth term (J2i=o "') ls zer0 ^ 
k = 2 or k = 3, and the notations 7 and 5 are slightly abusive : 7 stands for the function 

(u,..., u (k - 3 \x, y,z,v,..., v {e ' 1} ) ^ j( x , y, z, a(u, u {k - 3 \x, y,z,v,..., ^~ 1} )) , 

and the same for 5. With the same abuse of notations, (16-e) reads 

X- f - S ^ 0. (83) 

The equalities (g g^f-i) + g J_ 1} )u^ k ~ 2 ^ = ^w (fc ~ 2) = du{ 9 k -D u {k ~ 2) = are obvious in the 
original coordinates. Since the inverse of the change of coordinates (80) is given by 

= X (u, u( k -V ,x,y,z,v,..., = Y X (u, . . . , , x, y, z, v, . . . , v^), 

and E, X and X\ are given by (77), those equalities imply 

E X = X X = X lX = . (84) 

Then, from (75), (81) and (82), 

X 1 = {X 1 -8)^- + (X5-Ya)-?- , (85) 
ay oz 

[X, Xi] = (X 2 1 - 2X8 -Ya)— + (X 2 S - XYa - X t a) — . (86) 

oy oz 
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With these expressions of X and X\, the first relation in (78) implies : 

. (87) 



X-y-S X 2 -/ - 2X5 + Ya 
XS-Ya X 2 6 - XYa- X ia 



The definition of a implies Xz — a. In the original coordinates, this translates into 
the identity p xx = a(u, . . . ,u {k ~ 3, \x,p,p x ,v, . . ^v^' 1 **). Since Ep = Ep x = Ep x . x = (see 
(20)), applying E to both sides of this identity yields that a does not depend on its argument 
Also, if k > 3, differentiating both sides with respect to u^ x \ u^ -2 ) and u (k ~ 3 \ 
we obtain that the determinant (37) is zero if and only if a does not depend on its argument 
M ( fc -3) t q gum U p^ un( j er the assumptions of the lemma, 

a depends on u, . . . , u^ k ~ 4 \ x, y, z, v, . . . , v^~ 2 ^ only (88) 

with the convention that the first list is empty if k = 2 or k = 3. Now define two integers 
m < k — 4 and n < I — 2 as the smallest such that a depends on u, . . . , u 1 - 171 ^ , x, y , v, . . . , v^ n \ 
with the convention that m < if k — 2, k — 3, or a depends on none of the variables 
u, . . . , i/ fe_4 ), and n < if a depends on none of the variables v, . . . , v^ 2 \ We have 

to > => a u( _m) ^ , n > => a v (n) ^ . (89) 

Since to is no larger that k — 4, \ does not appear in the expression of Ya : 

m n 

Ya = ja y + 5a z + J2 u(i+1)a u(*) + ^ (90) 
i=a i=o 

where the first (or second) sum is empty if to (or n) is negative. 

In the left-hand side of (87), all the terms depend only on u, . . . , u^ m \ x, y, z, v, . . . , , 
except Ya, XYa and X\a that depend on u^ m+1 ^ if m > or on if n > (see 

above); the determinant is a polynomial of degree two with respect to u( m+1 ) and 
with coefficients depending on u, . . . , u^ m \ x, y, z, v, . . . , i/™' only, and the term of degree 
two, coming from (Ya) 2 , is 

Hence (87) implies a u ( m ) — a„(n) = and, from (89), negativity of to and n are negative. 
By definition of these integers, this implies that a depends on (x,y, z) only: in the original 
coordinates, one has p xx — a(x,p,p x ). □ 

Before proving Lemma 3.12, we need to extract more information from the previous proof : 

Lemma B.l. Assume, as in Lemma 3.11, that p is a solution of S^'t satisfying (38), but 
assume also that i > k > 3 and the determinant ( 37) is nonzero. Then [X, E2] = [X, E3] = 0. 

Proof. Starting as in the proof of Lemma 3.11, one does not obtain (88) but, since (37) is 
nonzero, 

a depends on u, . . . , u( k ~ 4 \ x, y, z, v, . . . , i/ £ ~ 2 ' and a u (k-a) ^ . (91) 

Since Ep = Ep x = 0, one has E = d/dv^" 1 ^ in these coordinates. The first equation in (84) 
then reads Xv 1 - 1 - 1 ) = 0) an d (75) and (82) yield 

d d 
E2 =d^T)> [X,E 2 } = -a v «- 2) — . 

Since [X, E 2 ] = ^X x (see (79)), relations (85) and (83) imply that a v (e-2), n, and the bracket 
[X, E2] are zero, and prove the first part of the lemma. Let us turn to [X, E3] : from (75) 
and (82), one gets, since E 2 and X commute, and X\ = 0, 

S 3=X,(-,^|z3)+^3), [X,E 3 ] = -(E 3 a)§-. (92) 
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In order to prove that E^a = 0, let us examine equation (87). For short, we use the 
symbol O to denote any function that depends on u, . . . , u^ k ~ 3 \ x, y, z, v, . . . , v ( £ ~ 3 ) only. 
For instance, — S = O, and all terms in the determinant are of this nature, except the 
following three : 

Ya — x a «(M+»' H ' a t(«- ! )+^i 
XYa = x x ot u (k-3) + v {t ~ 2) Xa vi i-3) +0, 

X x a = -a 2 ^a„( t -3)+" (<_2) a#-3))+0 

(we used X\ — 0). Setting £ = x a u( k - 3 ~> + v^~ 2 ^ a v ^-3) , one has 

I( = ^ ( + bs N) with b = ^ M) -^-3,^, (93) 
a u ( k -3) a u ( k -3) 

and equation (87) reads 

C 2 + OQ~{X 1 -8)hv^-^ +0 = 0. (94) 
Differentiating with respect to X and using (93) yields 

„ Xa u ( k -3) 



. 



a u {k-3) 

Then, eliminating £ between these two polynomials yields the resultant 

1 O -(Xj-S)bv( £ -V +0 

1 O -(X 7 -(5)W^- 2 ) +o 

2 X "" (fc - 3) 2b^- 2 )+0 0/- 2 >+0 

2 X "" (fc - 3) 2bv^+0 OvV-Q+O 

This is a polynomial of degree at most three with respect to v^ e ~ 2 ^ , the coefficient of (v^^ ) 3 
being — 4b 3 (X7— 5). Hence b = and, from (94), £ does not depend on y( e ~ 2 \ This implies 
E 3 a — because, from (92) and the definition of £, one has Cv('- 2 ) = E^a. □ 

Proof of Lemma 3.12. The independent variables in £33 are u 7 u,u,x,v,v,v. Since the de- 
terminant (37) is nonzero, one defines local coordinates (x, y, z, w, v, v, v) by 

y = p(u,u,u,x,v,v,v), z = p x (u,ii,u,x,v,v,v), w — p xx {u,u,u,x,v,v,v). (95) 

In these coordinates, X and Y, defined in (74), have the following expressions, with tp and 
a some functions to be studied further : 

d d d d 
X = — + z —+w—+a—, 96 
ox ay oz ow 

d . d , d . d .. d 
oy oz Ow Ov Ov 
Then, using, for short, the following notation T : 

T = Xj-6 + , (98) 

one has 

Xl =T^ + (XS-^§- z + (X^-Ya)-^, (99) 

f) f) f) 

[X, Xi] = (XT — XS + ip) — — h (X 2 S - 2A> + Ya) — + (X 2 ^ - XYa - X x a) 7^. (100) 

oy ' Oz v ' ow 

Also, 

E=±, E 2 = [E U Y]=^ + ^, [X,E 2 ] = I**- + (X^ - E 2 a) °- 
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but, from Lemma B.l, one has [X, E 2 ] = 0, hence ip$ = 0, E 2 = d/dv and — 0. Then 
E 3 = [JU] = + [X,Ek] = ^ + (X^ - E 3 a) A , 

but, from Lemma B.l, one has [X, E 3 ] = 0, hence ipi, = 0, £3 = d/dv and a v — 0. To sum 
up, 

E =7F> E 2 = -^, E 3 = ^-, (101) 
dv dv dv 

a depends at most on (x, y, z, w) only and ip on (x, y, z, w, v). 

Notation: until the end of this proof, O stands for any function of x,y,z,w only. For 
instance, a = O, 7 = O, 5 = O, T = O, XT = O , XS = O and X 2 S = O. 

r xs-tp 

XT — XS + ip X 2 5 - 2Xijj + Ya 
Xip = ^-iP 2 + OtP + (D . (102) 



From (78), (99) and (100), one has 



. Hence 



2T 

We now write the expression (99) of X\ as 

X x = X° + iPX\ + ^X\ (103) 

Note that X®, X\ and X\ are vector fields in the variables x, y, z, w only. Now define : 

p --«-fS-5 + H)s' (1 » 5) 

v _ x?-^? _ r » +0 » + ( -£)«., (106, 

so that 

X 1 = V-ipU (107) 

and, from (97) and (103) one deduces the following expression of X2 = : 

X 2 . mV+(X ^)l + ^l^(^l i + 0§- z+ 0l) + ^l + X° (108) 

where X\ and X® are two vector fields in the variables x, y, z, w only. 

This formula and (101) imply [£2, X2] = (Yip}^ U — ip v U; hence, from the last relation 
in (79), either ip v is identically zero or U is a linear combination of X\ and X 2 . We assume, 
until the end of the proof, that U is a linear combination of X\ and X 2 . This implies, using 
(107), that X 2 and X\ are linear combinations of U and V; hence U, V is another basis for 
X\,X 2 . Also, from (78) [U,V] must be a linear combination of U and V. From (105) and 
(106), 

where W 1 and W° are two vector fields in the variables x,y,z,w only, and, finally, with Z 1 
and Z° two other vector fields in the variables x, y, z, w only, one has, from (108) 

This vector field is also a linear combination of U and V. Computing the determinant in the 
basis d/dy, d/dz, d/dw, one has, using (105) and (106), 

det (U, V, X 2 - (Yi>)U - T[U, V] ) = ^ A + Oip 3 + Oip 2 + Oip + O = 0. 

It is assumed from the definition of that the partial derivative of 7 with respect to its 
fourth argument is nonzero; hence j w ^ and the above polynomial of degree 4 with respect 
to ip is nontrivial; its coefficients depend on x, y, z, w only, hence ip cannot depend on v. 
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We have proved that, in any case, both a and ip depend on x, y, z, w only, and this yields 
the desired identities in the lemma. □ 

Appendix C 

Lemma C.l. Let p be a smooth function of u, . . . , u^ 1 ^, x, v, . . . ,v^ _1 \ r a smooth func- 
tion of u, . . . , u^^^jX, v, . . . , with r v (t) ^ 0, and f a smooth function of four variables 
such that 

fc-2 l-l 

J2 u(l+1) Pu(*> + rp u(k -D + ^V i+1) p„w = f( x >P>P^Pxx) ( 109 ) 

where, by convention, rp u (k-i) is zero if k — and the first (resp. last) sum is zero if k < 1 
(resp. i — 0). Then either p depends on x only or 

k>l, £>1, p u (*-i)^0, p vlt -i)?0. (110) 

Proof. Let to < k — 1 and n < I — 1 be the smallest integers such that p depends on 
u, . . . , u^ m \ x, v, . . . , t/ n ); if p depends on none of the variables u, . . . , u^ 1 ^ (orv,..., u^ -1 )), 
take m < (or n < 0). Then p„(m> 7^ if m > and p^i^) 7^ if n > 0. 

The lemma states that either m < and n<0orfc>l,£>l and (m, n) = (k — l,£— 1). 
This is indeed true : 

- if m = fc — 1 and k > 1 then n = £ — 1 and £ > 1 because if not, differentiating both sides in 
(109) with respect to would yield r v (e)p u (k-i) — 0, but the lemma assumes that r v w ^= 0, 

- if m < fc-1 or m = 0, (109) becomes : M (l+1) p„(«) +E"=o ^ l+1) p v («) = f(x,p,p x ,p xx ); 
if to > 0, differentiating with respect to u( m+1 ) yields = and if n > 0, differentiating 
with respect to u( m+1 ) yields = 0; hence to and n must both be negative. □ 
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